Simplified Hydrostatic Carbon Burning in White Dwarf Interiors. 

Francisco Forster 1,3 Pierre Lesaffre 2 and Philipp Podsiadlowski 3 

Received ; accepted 

O ' 
(N 

Oh' 

' Not to appear in Nonlearned J., 45. 

C/3 ' 
00 

p? 

on 

6 
& 

> 

(N 
VD 

ON 

o 
o 



X 



1 Departamento de Astronomia, Universidad de Chile, Casilla 36-D, Santiago, Chile 
2 Laboratoire de Radioastronomie, 24 rue Lhomond, 75231 PARIS Cedex 05, France 
3 University of Oxford, Department of Physics, Oxford, 0X1 3RH, UK 



- 2 - 
ABSTRACT 

We introduce two simplified nuclear networks that can be used in hydrostatic 
carbon burning reactions occurring in white dwarf interiors. They model the 
relevant nuclear reactions in carbon-oxygen white dwarfs (CO WDs) approaching 
ignition in Type la supernova (SN la) progenitors, including the effects of the 
main e~-captures and /3-decays that drive the convective Urea process. They 
are based on studies of a detailed nuclear network compiled by the authors and 
are defined by approximate sets of differential equations whose derivations are 
included in the text. The first network, Nl, provides a good first order estimation 
of the distribution of ashes and it also provides a simple picture of the main 
reactions occuring during this phase of evolution. The second network, N2, is a 
more refined version of Nl and can reproduce the evolution of the main physical 
properties of the full network to the 5% level. We compare the evolution of 
the mole fraction of the relevant nuclei, the neutron excess, the photon energy 
generation and the neutrino losses between both simplified networks and the 
detailed reaction network in a fixed temperature and density parcel of gas. 

Subject headings: supernova, white dwarfs, nuclear reactions 
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INTRODUCTION 



Type la supernovae (SNe la) are the thermonuclear explosions of white dwarf stars. 
They are observable end-points of stellar evolution, they shape the energy and chemistry 
evolution of galaxies and have been successfully used as distance i ndicators up to redshifts 



of ~ 1.7 thanks to an empirical decline rate - luminosity relation ( jPskovskii 



1977 



1998 



1993) . This relationship led to the discovery of the acceleration of the Universe ( jRiess et al. 



Phillips 



Perlmutter et al. 



19991 ) 



The exact nature of SN la progenitors is still debated (IHillebrandt &: Niemeyer 

20001 ). but most models involve an accreting c arbon-oxygen white dwarf (CO WDs ) 



19841 ) 



with a mass close to the Chandrasekhar mass (jNomoto. Thielemann, fc Yokoi 
This CO WD would accrete mass from either a slightly evolved main sequence 
companion (CO WD + MS), a more evolved red giant star (CO WD + RG), or a 



Hach 


isu. Kato. & Nomoto 


„ , , „„ 
1996; 


Li &; van den 


ieuvel 


1997: 


Hachisu. Kato. & Nomoto 




1999: 


Hachisu. Kato. Nome 


)to. & 


Jmeda 


1999: 


Lander. Deu 


;schmann. Wellstein. HoefTich 


2000; 


Han & Podsiadlowski 


2004: 


Meng. Chen. & Han 


2009 


), or as a result of the merger 



of two degenerate stars with a combined mass above t 



re Chandrasekhar limit, in the 



double degenerate scenario (DD, see e.g. 



Webbink 



1984 : 



Iben fc Tutukov 



1984 ). It has 



been suggested that only CO W Ds accreting within a na r row range lead to successful 



therm onuclear explosions (see e.g. 



Nomoto fc Kondo 



1991 



Nomoto. Saio. Kato. fc Hachisu 



20071 ). unless additional physical p rocesses are included in t 



2007 



-dimensional e 


fects 


Pakmor et al. 


2010) 



Dommguez et al. 



2006 



re models, e.g. rotation and other 



Yoon. Podsiadlowski. &: Rosswog 



The diversity of SN la ejecta and the origin of the decline rate-absolute magnitude 
relation are being understood only recently thanks to new observational techniques 
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and energy transfer codes (see e.g. 


M 


azzali 


. Ropke. Benetti. & t 


illebrandt 


2007 




Wooslev. Kasen. Blinnikov. & Sorokina 


2007 


Kromer & Sim 


2009 


). These efforts 


have 



been accompanied by new developments in the physics of the explosion which have shown 
that even if pure d eflagration models can reproduce many SN la spectra a nd light curves 



(Ro pke et a 



Ropke et al. 



2007J), in some cases a delayed detonation may be necessary (IKhokhlov 



1991 



20071 ) . The physics of b oth the transition to detonation and ignition of th e 



Ropke 



2007 



wave remain uncertain ('. 


"apichino. Briiggen. Hillebrandt. & Niemever 


2006 


Iapichino & Lesaffre 


20K 


])• 



A related problem, rarely addressed in the literature, is how to connect theoretical 
models with obser ved systematic diffe r ences between SNe la occurring in different stellar 



environments (e.g. 



Hamuy et al. 



1996 



Sullivan et al. 



20061 ). These should be related to 



the pre-supernova evolution and no t to line-of-sight effects or othe r processes occuring 



randomly during the explosion (e.g. 



Kasen. Ropke. &: Wooslev 



2009 



Maeda et al. 



201fl). 



'hus, presupernova evolution must play a significant role in the diversity of SN la explosions 
(jLesaffre et all 12006 ). 



1.1. Presupernova Evolution: from Cooling to Ignition 



Here we will assume that the progenitors of SNe la originate in the SD scenario; 
but note that, even in the standard DD scenario, the core would evolve in a very 
similar way in the immediate pre-expl osion phase (i.e. during the last ~ 10 3 yr; see 



Yoon. Podsiadlowski. fc Rosswog 



20071 1 . 



Before a SN la progenit or becomes unbound by the explosion it undergoes sever al 



distinct phases of evolution ( INomoto. Thielemann. &: Yokoi 



1984 



Lesaffre et al. 



20061 ) 



First, the progenitor white dwarf cools down at almost constant density after its birth, 
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for a period of hundreds or thousands of Myr, depending on the particular formation 
scenario (the cooling phase). Then, it accretes matter for a period of the order of one 
Myr, making the degenerate star shrink to keep the hydrostatic equilibrium and its core 
compress adiabatically, the accretion phase. Adiabatic compression at the center and 
heat diffusion from the hot accreting envelope can make the central temperature increase 
under degenerate conditions, triggering hydrostatic carbon burning and a new source of 
energy generation. The energy input from hydrostatic carbon burning will force the star to 
transport the excess energy at its center convectively, the simmering phase. During this 
phase the convective core can grow to engulf most of the star. If the central density is high 
enough, e~-captures and /3-decays can become important. In the presence of a convective 
core, that process has been called the convective Urea process. 

At some point energy deposition will dominate over the energy losses, making the 
star's central temperature increase at almost constant density, the thermonuclear flash. 
Finally, when the temperature is high enough, one or more ignition spots will give rise to 
nuclear flames that will propagate in the highly convective medium, the thermonuclear 
runaway, causing a deflagration wave to sweep over the star, sometimes transitioning into a 
detonation wave. The deflagration and detonation waves will generate most of the kinetic 
energy and ashes in the ejecta in a few seconds, including radioactive matter which will 
later power the light curve of the supernova. Since this last phase will occur at very high 
temperatures, the characteristic burning time-scales will be much shorter than the typical 
weak interaction time-scales and the neutron excess of the ejecta will not differ from that 
of the presupernova star, except for the star's innermost regions where weak interaction 
time-scales are shorter. Most of the WD's neutron excess changes will occur before ignition. 
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1.2. Nuclear physics and the convective URCA process 



One of the main obstacles that remain to be solved in order to obtain self-consistent 
pre-supernova models up to ignition is the convective Urea process, which was mentioned 
in the previous section. The following factors make this phase of evolution difficult to solve: 
(1) the appearance of a fast-growing convective core with a very steep luminosity and 
composition gradient at its outer edge, the so-called C-flash; (2) a high-density medium 
with e~-capture and /3-decay time-scales similar to the convective time-scales, which have 
an uncertain effect over the energy budget of the star; (3) a steep density gradient which 
changes rapidly when the WD approaches the Chandrasekhar mass; (4) a steep composition 
gradient of the Urea pair 23 Na- 23 Ne at the threshold density for electron captures, which 
moves inwards as the central density increases and (5) an increasingly complex set of 
nuclear reactions as carbon burns and pollutes the WD with its ashes, with a time-scale 
that can be similar to e~-capture and /3-decay time-scales 

It is not clear whether e~-captures and /3-decays of Urea matter around a threshold 
density in a convective medium have a cooling or heating effect over time. Many stud ies 



over the years have reached different conclusions regarding this point. iPaczyhski 



( 119721 ) first 



suggested that Urea processes have a stabilizing effect over carbon burning , lead ing to the 



formation of a neutron star instead of a thermonuclear explosion. 



Bruennl ( 119731 ) realized 



that e -captures can cause heating by creating holes in the Fermi sea, which dom inate 



Couch fc Arnett 



(119741 ) found 



over the neutrino losses for the most important Urea pairs, 
that significant work must b e done to mantai n cony ection when Urea processes occur, 
with a net cooling effect, but iRegev fc Shavivl (|1975l) point e d out that convection cannot 



develop fast enough to prevent heating. iBarkat fc Wheeler! ( 119901 ) summarized the factors 



controlling the Urea process, but it was later shown that t 



should have been included in their analysis ( iMochkovitch 



le ro l e of the kinetic en ergy flux 



1996 



Stein et al. 



19991 ). More 



recently, iLesaffre. Podsiadlowski. fe Toutl ( 120051 ) showed that the heating effect of the Urea 



process depends on the state of mixing of the convective core, that the convective velocities 
are reduced by Urea processes and that time-dependent computations with a full nuclear 
network are needed to understand the effect of Urea processes on the ignition conditions of 
SNe la. 

In this study we have focused on how to accurately treat the increasingly complex 
nuclear reactions as the WD is polluted with ashes. We do not attempt to answer whether 
Urea processes cause cooling or heating, but provide a tool for better evaluating the 
competing heating and cooling mechanisms in a presupernova WD approaching ignition. 

In what follows we will introduce two approximations that use a limited number of 
nuclei to accurately describe the evolution of the main species that result from the burning 
of 12 C, the energy deposition rate, the energy losses via neutrinos and the neutron excess. 
We will show when these approximations hold, and their potential applications, but first 
we will introduce a detailed nuclear network which will be used for comparison with the 
simplified networks. 



2. THE FULL NUCLEAR NETWORK 

The nuclear reactions within a CO WD approaching explosion are characterized by 
the burning of 12 C nuclei at high densities (> 10 7 g cm -3 ) and high temperatures (> 10 s 
K) in an environment rich in 12 C and 16 O nuclei and relatively devoid of free protons, 
a-particles or neutrons. The dominant reactions are 12 C( 12 C, p) 23 Na, Q = 4.6 MeV, and 
12 C( 12 C, 4 He) 20 Ne, Q = 2.2 MeV, both occurring at similar rates (see Figured]). The 
nuclear network increases in complexity as the 12 C-burning pollutes the WD with its ashes, 
mainly 20 Ne and 23 Na, but also with protons and a-particles, which will be processed into 
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Fig. 1. — Ratio between the two main carbon burning reactions (thermally averaged cross- 
sections) plotted against the temperature in 10 9 K. 

additional 16 and 13 C nuclei, as will be shown later. 

We have integrated a detailed reaction network at fixed temperature and density 
compiled by the authors. The system of differential equation s defining the network is 



Bader fc Deuflhardl ( 119831 ). 



solved for using the semi-implicit extrapolation method from 
We include species that are part of the reactions known to be most important, or those 
close to them in a table of nuclides. We includ e all reaction rates betw een species of our 



network that are eit her in the Reaclib library ( [Thielemann et al 



interaction ta 



Dies 



Zegers et al 



^y 



Oda et al 



20061), or in the weak 



(11994]), as well as improved 13 N(e ,u e ) l3 C rates from 



( 120081 ) . Screening corrections were also included under the simplification 
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of iGraboske. Dewitt. Grossman. &: Cooperl (119731 ). Recent formalisms that treat carbon 
burning and screening corrections in alternative ways were not implem ented for this work 



^see e.g. 



Gasques et al. 



2005 



Yakovlev et al 



2006 



Spillane et al 



20071 ). 



We have assumed that the initial C/O ratio is always given by the nuclide mass 
fractions X( 12 C) = 0.3 and X( 16 0) = . 7, in order to compare with the work by 
Chamulak. Brown. Timmes. fc Dupczaki ( 120081 ). but reasonable variations of the initial 
compositions do not affect the validity of our approximations. Figure [2] show the main flows 
in our nuclear network at a temperature of 3 x 10 8 K and densities of 1 x 10 9 g cm -3 (top) 
and 3 x 10 9 g cm -3 (bottom), starting from a pure carbon-oxygen mixture, when 20% of 
the original carbon nuclei have been burnt. The left hand side color-coding bar indicates 
the mole fraction scale of the individual species, whereas the right hand side color-coding 
bar indicates the mole fraction flow scale of the individual reactions. Only reactions that 
are bigger than one thousandth of the biggest flow are plotted. 

The main differences between the low density (top) and high density (bottom) flow 
diagrams are due to e~-captures in the 23 Na(e~, z/ e ) 23 Ne reaction, which only occurs at 
densities higher than ~1.7xl0 9 g cm -3 . Below this density the e~-captures are almost 
exclusively due to the 13 N(e~, z/ e ) 13 C reaction. 

The advantage of this detailed network is that it can accurately track the evolution 
of nuclei with a wide range of characteristic time-scales. For example, it is capable of 
accurately following the mole fraction of a-particles, protons, neutrons or fast electron- 
capturing 13 N nuclei, as well as long characteristic time-scale nuclei like 12 C, or intermediate 
time-scale nuclei like 23 Na at high densities. 

Since we study this reaction network in convective WD interiors, we must consider the 
relation between the relevant burning time-scales and the convective turn-over time-scale, 
which is determined by the diffusion time-scale and the temperature, density and 
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Fig. 2. — Main flows at a temperature of 4xl0 8 K and densities of lxlO 9 g cm 3 (top) and 
3xl0 9 g cm -3 (bottom) when 20% of the original carbon has been burnt. 
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composition gradients of the star. The convective time-scale will be the shortest time-scale 
that affects the evolution of the structure of the star that is not directly co nnected to the 
nucle ar physics, and can be as low as 50 sec before ignition (see Figure 6 



Lesaffre et al. 



200J). 



Ideally, we would like to build an approximated nuclear network where we remove 
time-scales much shorter than the convective time-scale from the resulting set of differential 
equations, as long as a reasonable physical justification is provided. This would reduce the 
number of variables per zone while keeping accuracy and would make the inversion of the 
corresponding Jacobian in a Newton-Raphson integration scheme more stable. 

We have found two approximate solutions that achieve the former based on the 
deta iled nuclear network at fixed density and temperature described in t his section (see 



also 



Chamulak. Brown. Timmes. fc Dupczak 



2008 



Piro fc Bildsten 



20081 ). The first such 



approximation will be described in the following section. 
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3. THE SIMPLIFIED NETWORK: FIRST APPROXIMATION (Nl) 



To accurately follow the chemistry of the evolving WD progenitor towards explosion, 



the s imultaneous solution of the structure and chemistry of the star is required (jStancliffe 



20061 ). The chemistry equations must be able to reproduce the evolution of the main 
species and their effect on the star through changes in energy deposition, energy losses or 
the electron fraction. This can become computationally demanding if too many species are 
included a nd can worsen du ring the simmering phase of evolution due to known convergence 



problems (jlben 



1978 



19821 ). even with approximate theories of convection. 



We have found a first order approximation to the full nuclear network which will be 
described in this section. We will refer to this approximation as Nl. A more accurate version 
of this approximation, N2, will be described in Section [5j Both simplified networks use 
the fact that the dominant carbon burning reactions, 12 C( 12 C,p) 23 Na and 12 C( 12 C, a) 20 Ne, 
occur nearly at the same rate and the abundances of protons, a-particles, neutrons and 
13 N nuclei are approximately at equilibrium most of the time. This can be explained by the 
complementarity of the relevant reactions (see Figu re [3] and HI) , which was first noticed by 



Chamulak. Brown. Timmes. fc Dupczakl (120081 ) and 



Piro fc Bildstenl (120081 ) 



First, when in the carbon burning reaction 12 C( 12 C,p) 23 Na a proton is released, it will 
be captured quickly, preferentially in the 12 C(p, 7) 13 N reaction. Second, 13 N will quickly 
capture an electron in the 13 N(e~, t / e ) 13 C reaction, decreasing the pressure and lifting the 
degeneracy of the gas, and producing a 13 C nucleus. Third, in the carbon burning reaction 
12 C( 12 C, a) 20 Ne an a-particle will be released and quickly captured in the 13 C(a,n) 16 
reaction, thus consuming the 13 C nucleus just produced, but liberating a neutron. Fourth, 
the free neutron will be preferentially captured in the 12 C(n, 7) 13 C reaction, thus recovering 
the 13 C just consumed. Hence, the net effect is approximately the burning of six 12 C nuclei 




Fig. 3. — Main 12 C burning reactions and its ashes at low density. Six 12 C nuclei are 
consumed to produce four nuclei: 23 Na, 20 Ne, 16 and 13 C and one electron capture in the 
reaction 13 N(e~, z/ e ) 13 C. The figure is shown as a rotated tree, with parent nuclei connected 
to children nuclei in the direction indicated by the arrows. Double arrows correspond to 
electron captures or inverse beta decays. 

and the capture of one electron to be replaced by four nuclei: 20 Ne, 23 Na, 16 and 13 C, i.e.: 

6 12 C + e" -> 23 Na + 20 Ne + 16 + 13 C. (1) 

An additional e~-capture can occur via the 23 Na(e~, i/ e ) 23 Ne reaction if the density is above 
w 1.7 x 10 9 g cm -3 . In this case, two e - -captures can occur and the net effect is: 

6 12 C + 2 e~ -> 23 Ne + 20 Ne + 16 + 13 C, (2) 

which is schematically shown in Figure HJ If any of the flows above is broken the simplified 
network will fail and will need a different treatment. 

Now, we will assume that the only relevant reactions are: 



Fig. 4. — Same as Figure [31 but at high density. Six 12 C nuclei are consumed to produce 
four nuclei: 23 Ne, 20 Ne, 16 and 13 C and two e~-captures in the reactions 13 N(e~, ^ e ) 13 C 
and 23 Na(e-,i/ e ) 23 Ne. 

1) 12 C( 12 C,p) 23 Na, Q = 2.24MeV, 2) 12 C( 12 C, a) 20 Ne, Q = 4.62 MeV 

3) 12 C(p, 7 ) 13 N, Q = 1.94 MeV, 4) 13 N(e", i/ e ) 13 C, Q = 2.22 MeV 

5) 13 C(a,n) 16 0, Q = 2.22 MeV, 6) 12 C(n, 7 ) 13 C, Q = 4.95 MeV 

7) 23 Na(e-,z/ e ) 23 Ne, Q = -4.38 MeV, 8) 23 Ne(/3-) 23 Na, Q = 4.38 MeV. 

Note that only mass differences are used to compute the Q-values and not the energy 
of the electrons lost or gained in weak interactions and that the 13 N positron decay 
reaction must also be included at low densities. The change in electron density, which can 
be computed assuming charge conservation, can be used with the chemical potential of 
the electrons to compute the energy changes due to electron gains or losses. The Fermi 
energy of the electrons at high densities is app roximately 5.1 MeV [pV e /10 9 g cm" 3 ] 1 ' 3 



( IChamulak. Brown. Timmes. &: Dupczak 



2008|). 
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Since ignition can occur at temperatures as high as 8 x 10 8 K (see Figure 5 in 



Lesaffre et al. 



20061 ). it is necessary to consider whether any inverse reaction from the list 
above can be significant. We have also found that only the inverse of reaction 3), i.e. 
13 C(7,p) 12 C, needs to be considered, since its characteristic time-scale can be comparable 
or smaller than the characteristic time-scale of 13 N e~-captures, e.g. we have found that 
both time-scales are similar for a temperature of 8 x 10 8 K and a de nsity of 4 x 1 s g cm " 3 . 



Note that we use reaction rates for the N(e , u e ) 3 C reaction from 



Zegers et al 



ri2008h . 



Hence, we consider the following species: 12 C, 13 C, 13 N, 16 0, 20 Ne, 23 Na, 23 Ne, p, a and 
n. Defining \ = (av)i as the thermally averaged cross-section for the strong interactions 
1 to 6, or the rate of occurrence per particle per unit time per unit volume of the weak 
interactions 7 and 8, we can write the following differential equations: 



dY( 12 C)/dt = -Y 2 ( 12 C)pN A X 1 - Y 2 ( 12 C)pN A X 



dY( 13 C)/dt 
dY{ 13 N)/dt 
dY{ w O)/dt 
dY{ 20 Ne)/dt 
dY( 2 ^d)/dt 
dY( 23 Ne)/dt 
dY(p)/dt 
dY(a)/dt 
dY(n)/dt 



Y( 12 C)Y(p)pN A X 3 - Y( 12 C)Y(n)pN A X 6 + Y( 13 N)^\ 
Y( 12 C)Y(n)pN A X 6 - Y( 13 C)Y(a)pN A \ 5 + F( 13 N)A 4 , 
Y{ 12 C)Y(p)pN A X 3 - F( 13 N)A 4 - F( 13 N)A inv 



3 > 



YC 3 C)Y(a)pN A \ 5 , 
F 2 ( 12 C) 



2 

F 2 ( 12 C) 



pN A X 2 



pN A Xi - F( 23 Na)A 7 + F( 23 Ne)A 



F( 23 Na)A 7 - F( 23 Ne)A 

r 2 ( 12 c) 



2 

r 2 ( 12 c) 



pN A X, - YC 2 C)Y(p)pN A X 3 + F( N)A 



3 ' 



pN A X 2 - Y( Li C)Y(a)pN A X 5 



-Y( L2 C)Y(n)pN A X 6 + Y( Li C)Y(a)pN A X 



(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 

(11) 
(12) 



where X 3 nv is the thermally-averaged cross-section of the inverse reaction N(j,p) C. 
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Numerical experiments with a detailed network show that the mole fractions of protons, 
a-particles, neutrons and 13 N nuclei will be many orders of magnitude smaller than those 
of 12 C, 16 0, 20 Ne and, depending on the density and temperature, 13 C and 23 Na or 23 Ne. 
This means that their evolution is very fast compared to that of the rest of the species, 
even when very small quantities of 12 C have been burnt. Hence, we assume that the l.h.s. 
of equations (J5J), (1101) . (ITT]) and (lT2j) is much smaller than each individual term in the r.h.s. 
of their respective equations. Neglecting these time derivatives we can write the following 
equilibrium mole fractions for the nuclei p, a, n and 13 N, which will hereafter be referred 
to as the trace nuclei: 

m = n ,Q _^L /mv , y (a) = ^ 

Y(n) = Y( a C) A ?(13 N) = Y'Q'O ^1, (13) 

where / inv , defined as / inv = 1 + A^/A/i, indicates the relative strength of the inverse 
reaction 13 N(7,p) 12 C with respect to 13 N e~-captures. This factor is independent of the 
composition and we will normally have j mw 1, except if we approach ignition at relatively 
low densities, where 13 N e~-captures cannot compete with the inverse reaction 13 N(7,p) 12 C, 
e.g. in off-center ignition. 

If small quantities of 12 C are burnt, the mole fractions of the trace nuclei will reach the 
equilibrium values in equations ( 1T3|) . The typical time-scales for the equilibrium values to 
be reached from either lower or higher abundances can be found dividing equations ( Tl3l) 
by the positive terms in the r.h.s. of equations (fTOj) . ( TTTj) . (fl2|) and ([5]), or an arbitrary 
higher-than-equilibrium mole fraction by the negative terms in the latter equations. Both 
calculations give the same time-scales, assuming f inv = 1, namely: 

rip) = [Y C 2 C)pN A \ 3 ] -\ r{a) = [Y ( 13 C)piV A A 5 ] ~' 

r{n) = [Y ( 12 C)piV A A 6 ] ~\ r( 13 N) = Aj 1 . (14) 
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In contrast, the typical time-scale for 12 C burning will be 



r( 12 C) « [F( 12 C)piV A (A 1 + A 2 )] 



-i 



(15) 



and the typical time-scales for 23 Na electron-captures or 23 Ne /3-decays will be Af and 
A 8 1 , respectively. We compute these time-scales at the temperatures and densities relevant 
for white dwarf interiors (see Table [TJ and find that the following relations will be normally 
satisfied: 



Given that the convective time-scales found in pre-ignition white dwarfs will be 
normally bigger than the biggest trace nuclei time-scale shown in Table HJ it can be 
assumed that the trace nuclei will be in equilibrium even within moving convective eddies. 
The trace nuclei equilibrium mole fractions will be analogous to other equilibrium state 
variables like the temperature, which is defined by the assumption of local thermodynamic 
equilibrium inside every point of the star. In fact, in every stellar evolution code, as far as 
the authors are aware, energy diffusion during convection is computed assuming that local 
thermodynamic equilibrium is achieved in time-scales much shorter than the convective 
turn-over time-scale. 

Thus, we can assume that the trace nuclei equilibrium mole fractions are given by 
equations f fl3|) in every point of the star. Strictly speaking, the reactions needed to 
reach the equilibrium values will break the assumption of energy conservation under this 
approximation. However, since the trace nuclei equilibrium mole fractions will be very small 
compared to the 12 C mole fraction changes, this effect will be negligible during 12 C burning. 

It is worth noticing that the trace nuclei equilibrium mole fractions will change with 
the t( 12 C) time-scale. For example, when 13 C nuclei are synthesized by the burning of 
12 C, the equilibrium abundance of a-particles will decrease significantly as can be seen in 




(16) 
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Table 1: Characteristic time-scales for a typical CO WD interior composition when the 12 C 
burning time-scale is less than 10 10 yr. We assumed the lowest 13 C abundances found in 
simulations that start with solar abundance in order to compute r a . In all cases the trace 
nuclei characteristic time-scales are much smaller than r( 12 C), A7 , Ag" 1 , or the convective 
turn-over ti me-scale [50 se c befor e ignition at a density of 1.8 x 10 9 g cm -3 in the worst case 



according to 



Lesaffre et al 



(120061 )]. which validates the use of the approximated network in 



both radiative and convective energy transport regions. 
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equations ffTB"]) . Assuming that the mole fraction of 12 C is constant for this purpose, the 
a-particle mole fraction would decrease as d\nY(a) = — <ilnF( 13 C). 

With the assumption of trace nuclei equilibrium we can obtain a simplified set of 
equations which do not include terms with the trace nuclei typical time-scales. Replacing 
equations (ITS]) into equations ([3]), @, (JZJ), (]H]) and (JHJ), we obtain: 

dY( 12 C)/dt = -3 y2( ^ C W (Ai + A 2 ), (17) 
Y 2 ( l2 C) 

dY( 13 C)/dt = — — -pAT A Ai (18) 
Y 2 ( 12 C) 

dY( 16 0)/dt = l 2 ) P N A X 2 , (19) 

dY( 20 Ne)/dt = ¥^lpN A \ 2 , (20) 

dY ( 23 Na) /eft = Y2 ^ C \ n a X 1 - F( 23 Na)A 7 + F( 23 Ne)A 8 , (21) 

dY{ 23 Ne)/dt = y( 23 Na)A 7 - F( 23 Ne)A 8 , (22) 

which constitutes the system of equations for a simplified nuclear network. 

A quick inspection of these equations shows that only 12 C, 23 Na and 23 Ne need to be 
tracked as primary species to accurately follow the chemistry changes in the star under this 
approximation, which will greatly simplify the computational cost of simultaneously solving 
the structure and chemistry of the star. Moreover, it can be seen that 12 C burns 50% faster 
than what one would naively obtain using only the two main carbo n burning reactions and 



igno ring the presence of small quantities of a s" 



and 



Chamulak. Brown. Timmes. fc Dupczakl (120081 ) 



ies 1 as first noticed in 



Piro & Bildstenl (I2008h 



We can also see that the evolution of the 13 C, 16 and 20 Ne nuclei will be initially 
faster than that of 12 C, but as significant amounts of 12 C are burnt their characteristic 
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time-scales will become comparable, i.e. assuming Ai ~ A2: 

r(-C) = 2 [rf^O^A,]- 1 „ 6 ^| r(»C), 

T(160) = 2 SI [>'( 12 C)P^A 2 ]-' « 6 T{ * C ), 
rfNe) = 2 [KCC^A.]"' „ 6 }g| r("C), 

The evolution of 23 Na and 23 Ne will be different and will depend on the density. Their 
characteristic time-scales will be: 

r(-Na) = min{2^p^[y(-C),iV, < A 1 ]- 1 , A 7 "\ £g| A 8 "'} (23) 

i(2iNe)=mm {Sj v ' v l (24) 

If the 23 Ne /3-decay time-scale were much shorter than the 23 Na e~-capture time-scale 
(A^ 1 <C Af 1 ) and the 12 C-burning time-scale [Ag 1 <C r( 12 C)], which corresponds to the low 
density limit and which we call hypothesis HI, 23 Ne nuclei produced by e~-captures from 
newly synthesized 23 Na would move to an equilibrium value in a time-scale Ag . This value 
would be obtained neglecting time derivatives in equation (122]) : 

f( 23 Ne) = y( 23 Na)^, (25) 
As 

and using this value in equation ( 12~T1) we obtain: 

Y 2 ( 12 C) 

dY( 23 Na)/dt = — — — -pN^Xi. (26) 

If the e~-capture time-scale were much smaller than the 23 Ne /3-decay time-scale 
(Af 1 <C Ag -1 ) and the 12 C-burning time-scale [Af 1 <C t( 12 C)], which corresponds to the 
high density limit and which we call hypothesis H2, 23 Na would act as trace nuclei and 
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with a time-scale A 7 1 would reach an equilibrium value, which could be obtained neglecting 
time-derivatives in equation f l2"Tj) : 

F(23 Na) = ^1 + F( 23 Ne) ^ (27) 

2 A7 A7 

Assuming this equilibrium value in equation ( )22|) we obtain: 

Y 2 ( l2 C) 

dY{ 23 Ne)/dt = \ > pN A X 1 (28) 

These alternative additional simplifications are useful for understanding the behavior 
of the simplified network in these extreme cases, but are not used in the simplification of 
this section. They are not valid if the 12 C-burning time-scale becomes smaller than both 
the 23 Na e~-capture and 23 Ne /3-decay time-scales, or if the latter time-scales are similar 
to each other. If this is the case, we must solve equations ( l2T|) and ( 1221) exactly under the 
simplified network. 



3.1. Dependence on the 12 C Burnt Mole Fraction 

We can also derive a set of differential equations that relate the increase in mole fraction 
of the secondary species with the amount of burnt carbon. First, we divide equations ( |T8l) . 
f lT9|) and flU by equation ^TTJ) to obtain: 

-0.15, (29) 
-0.19, (30) 
-0.19. (31) 



dY( 13 C) 


1 


dY( 12 C) 


3(1 + A 2 /Ax) 


dY( 16 0) 


1 


dY{ 12 C) 


3(1 + Ai/A 2 ) 


dY{ 20 Ne) 


1 


dY{ 12 C) 


3(1+Ai/A 2 ) 



For 23 Na or 23 Ne the result is not as simple, unless either of the conditions necessary 
for equations (126]) or (128]) are met. In those cases, we would get 



-22 - 



Under hypothesis HI (low density limit): 

dF( 23 Na) 



dY( 12 C) 3(1 + A 2 /Ai) 

Under hypothesis H2 (high density limit): 

dY{ 23 Ne) 1 



-0.15 (32) 



-0.15 (33) 



dY( 12 C) 3(1 + A 2 /Ai) 

3.2. Electron Mole Fraction 

Similarly, we can write equations for the evolution of the electron mole fraction Y e . 
Note that this quantity is closely related to the neutron excess, defined as r\ = Y^^iVu with 
7]\ = {n- x — pi)/ (n; + Pi) and where n; and p\ are the number of neutrons and protons of the 
respective species. The electron mole fraction and neutron excess are related by the formula 
T) = 1 — 2 Y e , which implies drj/dt = —2 dY e /dt. 

We note that when the 12 C-burning time-scale is much smaller than the e~-capture 
time-scale [r( 12 C) <C Xj 1 } and when 13 N is in equilibrium, the rate at which Y e changes will 
be the rate at which F( 12 C) changes due to the reaction 12 C( 12 C,p) 23 Na, i.e.: 

f = -rm pNAXl md _jp_ = _i_ M 0,5 . ,34) 

When the e~-capture time-scale becomes smaller than the 12 C burning time-scale, the 
e~-capture rate will be twice the former result, since 23 Na is produced in the same branch 
of the network as 13 N (see Fig. HJ). If both time-scales are comparable, the electron mole 
fraction can be obtained using that Y e = (1 — i])/2 and computing the neutron excess from 
the exact solution of equations ( TT71) to ( 1221) and the formula 



r) w F( 13 C) + F( 23 Na) + 3 F( 23 Ne) + Y(n) - Y(p) - F( 13 N). 



(35) 
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4. LIMITS OF THE SIMPLIFIED NETWORK Nl 

Now we will examine different ways the assumptions of the simplified network Nl can 
break down, and with this information build a second more refined simplified network, N2, 
which includes some of the leak reactions that will be discussed in what follows. 



4.1. Proton Leaks 

When the 23 Ne /3-decay time-scale is shorter than the 23 Na e~-capture time-scale, 
normally below the threshold density p t h ~ 1-7 g cm -3 , some of the protons that would 
be captured in the 12 C(p, 7) 13 N reaction can leak via the reactions 23 Na(p, a) 20 Ne and 
23 Na(p, 7) 24 Mg, which are defined as reactions 9) and 10), as the abundance of 23 Na 
increases. The ratio of their thermally averaged cross-sections is shown in Figure [5j 

Conversely, when the 23 Na e~-capture time-scale is shorter than the 23 Ne /3-decay 
time-scale, normally above the threshold density p t h ~ 1.7 g cm -3 , protons can leak via 
the reaction 23 Ne(p, n) 23 Na, which is defined as reaction 11), as the abundance of 23 Ne 
increases. These leak reactions can significantly change the distribution of ashes and energy 
input as 12 C is burnt, as well as the trace nuclei equilibrium mole fractions, as will be shown 
later. The abundances of 23 Na and 23 Ne at which the former reaction rates become equal 
to the 12 C(p, 7) 13 N reaction are shown in Figure El 

Let us define the cross-section for proton captures on either nuclei of the pair 
23 Na- 23 Ne, which we call Urea matter, as Au, and the cross-section for proton captures on 
12 C as Ac- The ratio between the rates in both reactions will be: 

_ y( 23 Na or 23 Ne) Au 
We can use the following approximate relation inferred from the simplified version of the 
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Comparison of proton captures onto Na 




0.2 0.4 0.6 0.8 

T9 

Fig. 5. — Cross-section ratio for proton captures onto 23 Na. The dominant proton-leak 
reaction will be 23 Na(p, a) 20 Ne for all the temperature range considered here. 



- 25 - 



Critical number ratio for proton captures 



U 



X 




0.2 



0.4 



0.6 



0.8 



T9 



Fig. 6. — Mole fraction ratio between Urea matter and 12 C when proton captures onto Urea 
matter become equal to proton captures onto 12 C, plotted against the temperature in units 
of 10 9 K. Screening corrections are taken into account, using a density of lxlO 9 g cm -3 for 
23 Na and 2xl0 9 g cm" 3 for 23 Ne, below and above the threshold for electron captures. 
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network: 



dY{ 



23 Na Qr 23 Ne ) 
dY( 12 C) 



0.15 



(37) 



and defining / as the burnt fraction of 12 C, we can write: 



r 



-0.15 AY( 



( 12 C) An 
) A c 




(38) 



F( 12 C) 



i.e. when r = 1, / « 6.7 Ac/Au, or 6.7x the ratio shown in Figure EJ Thus, when 9% or 
4% of the original 12 C is burnt, depending on whether protons leak on 23 Na or 23 Ne, proton 
leaks will become significant. From Figure [6] it can be noted that proton leaks will be 
stronger at temperatures of about 3 x 10 8 K. Above temperatures of 5 x 10 8 K the fraction 
of burnt 12 C for proton leaks to be important will change to approximately 27% and 14%, 
when either the 23 Ne /3-decay time-scale or the 23 Na electron capture time-scale is shorter, 
respectively. 

Although proton leaks start later when the 23 Ne /3-decay time-scale is shorter, their 
effect in this case is more difficult to model. Since less 13 N will be present due to the 
bypassing of one of the branches of the network, the amount of e~-captures is reduced and 
less a-capturing 13 C synthesised, opening secondary channels for a-captures. With more 
proton leaks, 12 C and 16 can become the main targets for a-captures, depending on the 
temperature as will be discussed later. 

When the 23 Na e~-capture time-scale is shorter, proton leaks will start sooner, which 
does not change significantly the main results of the simplified network because the end 
products will be the same after two e~-captures and the burning of six 12 C nuclei: 



The main difference will be that now both e -captures will be on 23 Na rather than on 13 N 
and 23 Na, changing the typical time-scale of the network in this leak branch. 



6 12 C + 2 e~ 23 Ne + 20 Ne + 16 + 13 C. 



(39) 
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The ratios between the relevant thermally averaged cross-sections for proton captures 
at different densities are shown in Figure [71 We summarize the proton-leak reactions in 
Figures El and 



+ 




0.2 



0.4 



0.6 



0.8 



T9 



Fig. 7. — Ratio between thermally averaged cross-sections relevant for proton captures at 
different densities [gr cm -3 ]. Solid lines correspond to (Ag + A10VA3 and dashed lines, to 
A13/A3 (see Section POT) . The threshold density for e~-captures onto 23 Na is close to 1.7xl0 9 
g cm -3 . The density dependence is due to screening corrections. 
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Fig. 8. — Same as Figure El but with proton leaks (dotted arrows). See discussion in the 
text for neutron and a-leaks. 

4.2. Neutron Leaks 

A similar analysis can be done to compute the critical mole fraction of secondary 
nuclei for neutron leaks to be important. In this case, the most important neutron capture 
reactions are 12 C(n, 7 ) 13 C, 20 Ne(n, 7 ) 21 Ne, 21 Ne(n, 7 ) 22 Ne and 23 Na(n, 7 ) 24 Na. The 20 Ne and 
21 Ne neutron capture reactions have very similar cross-sections, approximately seven times 
the 12 C neutron capture reaction. The 23 Na reaction has a cross-section approximately 11 
times bigger than the cross-section of the 12 C neutron capture reaction. 

The ratio of the mole fractions of 20 Ne and 23 Na with respect to 12 C can be related to 
the mole fraction change of 12 C using equations (13T1) and (1321) . at densities when the 23 Ne 
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Fig. 9. — Same as Figure HJ but with proton leaks and the resulting bypass for neutron 
captures (dotted arrows). See discussion in the text for neutron and a-leaks. 

/3-decay time-scale is shorter than the 23 Na e~-capture time-scale: 

« -0.34 AF( 12 C), (40) 

or when the 23 Na e~-capture is shorter than the 23 Ne /3-decay time-scale: 

F( 20 Ne) + F( 23 Na) « -0.19 AF( 12 C). (41) 

In the first case, using equation (jlO]) and using the added cross-sections of 20 Ne, 21 Ne and 
23 Na neutron captures, it can be shown that neutron leaks, i.e. significant neutron captures 
on species different than 12 C, occur when the fraction of burnt carbon is approximately 12% 
its original amount. In the second case, using equation (|41|) it can be shown that neutron 
leaks will occur when about 38% of the original carbon is burnt. 
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As first noticed by 



Chamulak et al 



(120071 ). neutron captures onto Fe will be 



negligible, since the 56 Fe(n, 7) 57 Fe reaction has a cross-section approximately 64 times 
bigger than that of neutron capture onto 12 C, but the mole fraction of 56 Fe is approximately 
1250 times smaller than that of 12 C at solar metallicity in the 12 C-rich environment of a 
WD. Thus, a 56 Fe abundance of more than 20 times the solar abundance would be required 
to compete with the reaction 12 C(n, 7) 13 C. 



4.3. a— particle Leaks 

Our detailed nuclear network shows that the most important reactions for a-captures 
are 13 C(a,n) 16 0, followed by the much weaker reactions 12 C(a,7) 16 and 16 0(a, 7) 20 Ne. 
In Figure [10] we show the ratios of the mole fractions of 12 C and 16 relative to the mole 
fraction of 13 C necessary for 13 C(a,n) 16 to be the dominant a-capture reaction. It can be 
seen that only a minimal amount of 13 C is needed for this to be the case: Y ( 13 C) ^3x 10 -8 , 
which is ~ 10 4 times lower than the solar metallicity value. 

Moreover, using arguments similar to those used for proton and neutron leaks, a 
fraction of only 4 x 10~ 6 of the original 12 C needs to be burnt to reproduce this abundance 
with zero metallicity, assuming that for every six 12 C nuclei one 13 C nucleus is produced. 
Thus, under the temperature range investigated here, one would conclude that these 
secondary reactions are negligible. 

However, if significant neutron leaks occur, 13 C can be significantly depleted and these 
reactions can become important. In fact, after 13 C is depleted as a consequence of neutron 
leaks and the a-particle over-production accompanying the 23 Na(p, a) 20 Ne proton leak 
reaction, a-leaks are necessary to reproduce the abundances of 16 0, 20 Ne and to limit the 
growth of a-particles and reproduce the 13 C depletion more accurately. 
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Critical number ratio for alpha captures 




0.2 0.4 0.6 0.8 

T9 



Fig. 10. — Number ratio of 12 C or 16 to 13 C when the rates of the a-capture reactions 
12 C(ci!,7) 16 or 16 0(o;, 7) 20 Ne equals the rate of a-captures in the reaction 13 C(o;, n) 16 0. 
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5. THE SIMPLIFIED NETWORK: SECOND APPROXIMATION (N2) 

We can now include the leak reactions discussed in the previous section and additional 
ones with the following indices: 



9) 


23 Na(p, a) 20 Ne, 


Q = 


2.38 MeV, 


10) 


23 Na(p, 7 ) 24 Mg, 


Q = 


11.69 MeV, 


11) 


23 Ne(p,n) 23 Na, 


Q = 


3.59 MeV, 


12) 


20 Ne(n, 7 ) 21 Ne, 


Q = 


6.76 MeV, 


13) 


23 Na(n, 7 ) 24 Na, 


Q = 


6.96 MeV, 


14) 


12 C(a, 7 ) 16 0, 


Q = 


7.16 MeV, 


15) 


16 O(a, 7 ) 20 Ne, 


Q = 


2.84 MeV, 


16) 


21 Ne(n, 7 ) 22 Ne, 


Q = 


10.36 MeV, 


17) 


13 C(p, 7 ) 14 N, 


Q = 


7.55 MeV, 


18) 


21 Ne(p, 7 ) 22 Na, 


Q = 


6.74 MeV, 


19) 


20 Ne(a, 7 ) 24 Mg, 


Q = 


9.32 MeV, 


20) 


23 Ne(a,n) 26 Mg, 


Q = 


5.41 MeV, 


21) 


23 Na(a,p) 26 Mg, 


Q = 


1.82 MeV, 


22) 


21 Ne(a,n) 26 Mg, 


Q = 


4.84 MeV, 



This list includes proton-leak reactions (9, 10, 11, 17 and 18), neutron-leak reactions 
(12, 13 and 16), and a-leak reactions (14, 15, 19, 20, 21 and 22). Proton-leaks are the first 
to be significant as 12 C ashes are produced, but with more ashes produced neutron-leaks 
become relevant too. When this happens, 13 C is rapidly depleted and a-leaks become 
important as well. 

As can be seen from the reactions in the list, it would be necessary to increase the 
number of species to account for all proton, neutron and a-leaks exactly. To do this we can 
either account for 24 Mg, 21 Ne, 24 Na, 22 Ne, 14 N, 22 Na and 26 Mg independently, or we can 
group some or all of them into an auxiliary variable. Since reactions 16) and 18) depend on 
the abundance of 21 Ne nuclei, we solve for 21 Ne independently and group 24 Mg, 24 Na, 22 Ne, 
14 N, 22 Na and 26 Mg into a dummy leak species. 

Thus, while keeping the number of independent variables small and ensuring mass 
conservation, we can introduce Y L as the mole fraction of a dummy nuclei that represents 
leaks on 24 Mg, 24 Na, 22 Ne, 14 N, 22 Na and 26 Mg. We assume that this dummy species has 



-33 - 



approximately the average mass number of the species it represents, defining it as: 

li^maxjo, ^|^ } (42) 

where A-^ is the mass number of the respective species and the sum is made over 12 C, 13 C, 
13 N, 16 0, 20 Ne, 21 Ne, 23 Ne and 23 Na. Since the leak nuclei 24 Mg, 24 Na, 22 Ne, 14 N and 
22 Na can subsequently capture neutrons with a cross-sections similar to that of 20 Ne, for 
simplicity we will assume that leak nuclei capture neutrons with the 20 Ne cross-section. 

If we now add new terms associated to the p, n and a-leak reactions discussed 
above into equations (j3J) to f|T2|) and assume a stationary solution for the trace nuclei 
p, a, n and 13 N, the modified equilibrium mole fractions can be described with the following 
equations (c.f. equations. [T3"|) : 

Y(p) = Y{p) / p , Y(a) = Y(a) f a , 

Y(n) = Y(n) f a f n , f( 13 N) = F( 13 N) / p , (43) 

where the auxiliary variable f p , f n and f a are defined as follows: 

f = L , F( 23 Na)(A 9 + Am) + r( 23 Ne)A n + F( 13 C)A 17 + F( 21 Ne)A 18 | 1 
/p -\ +/inv 1^( 12 C)A 3 J 1 ] 

y( 12 C)A 14 + F( 16 0)A 15 + r( 20 Ne)A 19 + F( 23 Ne)A 20 + F( 23 Na)A 21 + F( 21 Ne)A 22 
" F(i 3 C)A 5 

(45) 

/aSjr i/ 1 + 2 (46) 



L r 2 ( 12 c) a 2J 

= r F( 23 Ne) [y(p)A n + Y (a)A 20 ] + y( 21 Ne)l»A 22 

I + r( 13 C)F(a)A 5 
f [y( 20 Ne) + r L ]A 12 + y( 21 Ne)A 16 

I + F(i 2 C)A 6 



+ F( 23 Na)A 1 3^ 1 



(47) 

where we have neglected the contribution to f p from the leak reaction 23 Na(a,p) 26 Mg, which 
is generally small. An exact solution that takes into account the additional protons from this 
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or similar reactions can be obtained by multipliying the trace abundances Y(p) and F( 13 N) 
by the factor {1 + AY(a)} {1 - AY (a) BY (p)}' 1 , where A = 2F( 23 Na)A 2 i {F 2 ( 12 C)A 1 }- 1 
and B = 2^- 1 r( 23 Na)A 9 {F 2 ( 12 C)A 2 }- 1 . The factors f a and / n should be computed by 
using these modified values. 

These formulae can be used to approach the leak regime of the network and can be 
easily generalized to include more leak reactions involving protons, neutrons or a-particles. 
They predict a decrease in the number of free protons and 13 N nuclei at equilibrium and 
the fact that this decrease starts later when the 23 Ne /3-decay time-scale is shorter than 
the 23 Na e~-capture time-scale, at low densities. Also, they explain the contribution to 
a-particles from the proton leak reaction 23 Ne(p, a) 20 Ne, the a-leaks from 12 C and 16 
captures, the contribution to neutrons from the proton leak reaction 23 Ne(p, a) 20 Ne and the 
neutron depletion due to 20 Ne, 21 Ne and 23 Na neutron captures under our approximation, 
as well as a simple estimate of the neutron captures on 24 Mg, 24 Na, 22 Ne, 14 N and 22 Na. 

Now, we can rewrite the differential equations describing the evolution of the slowly 
varying nuclei 12 C, 13 C, 16 0, 20 Ne, 23 Na and 23 Ne, using the new trace nuclei equilibrium 



-35 - 



values and including the leak reactions discussed above, i.e. 



F(^C)piV A Y( u C)(\i + A 2 ) + + Y(n)A 6 + Y(a)A 14 , (48) 

I /inv J 

= Y( 13 N)A 4 + P iV A |r( 12 C)rHA 6 - Y( 13 C) [Y(a)X 5 + Y(p)X 17 ] J, (49) 

(/y 1 " 0j - piY A Y(a)(y( 13 C)A 5 + Y( 12 C)A 14 - Y( 16 0)A 15 1, (50) 



dt 

rfr( 20 Ne) Ar fy 2 ( 12 c) 



dt P nJ^^-X 2 + Y( 23 Na)Y(p)A 9 - Y( 20 Ne) [Y(n)A 12 + Y(a)X 19 ] 

+ Y( 16 0)Y(a)A 15 j, (51) 



dY( 2 Na) = -^( 23 Na)A 7 + Y( 23 Ne)A 8 + pN^^^X, + Y( 23 Ne)Y(p)An 

- Y( 23 Na) [Y(p)(X 9 + A 10 ) + Y(n)A 13 + Y(a)X 21 ] }, (52) 
dY( 23 Ne) 



Y ( 23 Na)A 7 - y ( 23 Ne)A 8 - Y ( 23 Ne)piV A [Y(p)X n + Y(a)X 20 ] , (53) 
P nJy ( 20 Ne)yHA 12 - y ( 21 Ne) [Y (p)A 18 + Y(n)A 16 + Y (a)A 22 ] ] . (54) 



dt 

rfY( 21 Ne) 



The energy generation rate can be obtained after a straightforward modification of 
the individual terms above, noting that the combined energy contribution of the inverse 
reaction 13 N(7,p) 12 C and the reaction 12 C(p, 7) 13 N can be obtained simply by multiplying 
the Q-value of reaction 3) to the associated term in the r.h.s. of equation (1481) . 



6. COMPARISON WITH THE DETAILED NETWORK 

In what follows we will compare the results of the full nuclear network introduced in 
Section [2] with those of the simplified nuclear networks Nl and N2 introduced in Sections [3] 
and |5j We will compare the evolution of the main and trace nuclei mole fractions, the time 
evolution and the energy release in the form of photons or neutrinos as a function of the 
fraction of burnt 12 C nuclei. 
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In most examples we have chosen a temperature of 4 x 10 8 K and a density of 3 x 10 9 
g cm -3 , which are typical of the conditions encountered during the thermonuclear runaway 
and before ignition in a pre-supernova CO WD. We chose a temperature and density at 
which the simplified solution errors would be representative of those encountered at other 
temperatures and densities after the trace nuclei have reached their equilibrium values. The 
time-scales for the trace nuclei to reach their equilibrium values can be computed from 
equations (fl4"|) . 

In Figures [TT1 and [T2l we show the ratio between the 12 C ashes and 12 C mole fraction 
changes vs the fraction of burnt 12 C, as well as the ratio between the neutron excess change 
and the 12 C mole fraction change vs the fraction of burnt 12 C. We compare the integration 
of the full network with that of the networks Nl and N2, respectively. The transformation 
was chosen because it makes the evolution of 13 C, 16 and 20 Ne look exactly flat under 
Nl, according to equations ([29]) . fl30l) and fl3"TT) . allowing an easier comparison between the 
integrations of Nl and N2. 

It can be seen that once the slowest evolving trace nuclei reach their equilibrium value 
and before about 1% of the 12 C has been burnt, both plots show a good agreement between 
the simplified and full networks for the main 12 C ashes. The exceptions are a-particles, 
16 and 20 Ne under Nl due to the initial a-leaks from reaction 15) and because we have 
assumed a pure CO mixture, with 13 C initially absent. Since in this example the 23 Na 
e~-capture time-scale is shorter than the 23 Ne /3-decay time-scale, both simplified solutions 
correctly predict a higher abundance for 23 Ne once both time-scales become comparable. 
Both solutions show a good match for the neutron excess changes too. 

However, once the fraction of burnt 12 C is above ~ 1%, the trace nuclei equilibrium 
abundances of Nl begin to differ from the exact solution, whereas N2 matches their values 
even when more than half of the original 12 C has been burnt. The inclusion of the leak 
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|AY( 12 C)/Y ( 12 C)| 

Fig. 11. — The ratio between the ashes and 12 C mole fraction changes vs the fraction of 
burnt 12 C in the full and Nl networks. Main ashes and the neutron excess change are shown 
as continuous (full network) or dot-dashed (Nl) lines and trace nuclei are shown as dashed 
(full network) or dotted (Nl) lines. See discussion in the text. 
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|AY( 12 C)/Y ( 12 C)| 

Fig. 12. — Same as Figure [HJ but using the network N2 as comparison. We also include 
the mole fraction changes of 21 Ne and the dummy leak particles introduced in equation (|42p . 
Note that the trace nuclei follow their equilibrium values precisely under this approximation. 
The abundance of leak nuclei is only well matched when more than pa 1% of the 12 C is burnt. 
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reactions in N2 provides a better match for the trace nuclei and as a consequence, the main 
12 C ashes as well, once the trace nuclei equilibrium mole fractions have been reached. Note 
also that the a-particle, 16 and 20 Ne mole fraction changes are well matched in N2 once 
trace nuclei equilibrium is reached. 

In N2 the trace nuclei follow precisely their equilibrium values as 12 C burns, even 
if these equilibrium values change dramatically. This is because their characteristic 
time-scales are much smaller than the characteristic 12 C burning time-scale. This supports 
the use of N2 in a varying temperature and density integration, which would be analogous 
to changing the equilibrium values of the trace nuclei as 12 C burns. The shortest time-scale 
for the environment's variables to change should be bigger than the biggest trace nuclei 
time-scale. 

In Figure [13] we show an enlarged section of Figure [121 If Nl were valid, the evolution 
of 16 and 20 Ne would look flat and approximately 0.19 in this space. The same would 
be the case for 13 C, although approximately 0.15 (see equations I2"9l [3TJ1 and l3"Tj) . This is 
only approximately true when the fraction of burnt 12 C is between 2% and 5%. Since we 
plot absolute values of mole fraction changes, Figure [13] does not show whether the 16 
mole fraction is decreased or increased. In fact, its mole fraction is originally depleted and, 
only after the fraction of burnt 12 C is about 10~ 4 , it is increased. The original depletion is 
caused by a-leaks in reaction 15), which increase the mole fraction of 20 Ne. 

The 23 Ne mole fraction is larger than the 23 Ne mole fraction only after the burnt 
fraction of 12 C is about 2 x 10~ 5 . Interestingly, this makes the neutron excess evolution 
change from above 0.3 to 0.6 in this transformation. This can be understood noticing that 
drj/dt = —2dY e /dt and that the value shown in equation is expected to double at high 
density. 

Note that 13 C is depleted as 21 Ne and leak nuclei increase their abundance. This 
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Fig. 13. — The ratio between the main ashes and 12 C mole fraction changes vs the fraction 
of burnt 12 C in the full (continuous) and N2 (dot-dashed) networks. Note that in Nl the 
evolution of the species 13 C, 16 and 20 Ne should be exactly flat in this plot according to 
equations (|29|) . (I30I) and (I3"IT) . Generally speaking, N2 reproduces the full network abun- 
dances to the 5% level when the trace nuclei have reached their equilibrium mole fractions. 
In numerical experiments, N2 loses its accuracy when both f p and f a (see equations HH and 
I4T1) are lower than 0.5, which is due to the observed 13 C depletion. 
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is because neutron leaks shortcut 13 C neutron captures, as discussed before. When this 
happens, f a decreases dramatically (see equation SB}, and the accuracy of the approximation 
is lost to about 20% level. Although at the beginning of the integration f a is also small, 
N2 matches the full network with great accuracy. To distinguish between the cases when 
the accuracy of N2 is very good and only a small fraction of 12 C has been burnt from the 
late loss of accuracy with a big fraction of 12 C burnt, we define the criterion for N2 to be 
considered an accurate representation of the full network as 

f p > 0.5 and f a > 0.5, (55) 

which should be used in a varying temperature and density integration, such as in real 
stellar evolution models. 

In Figure dH we show a similar integration as in Figure [T3] but at a density of 10 9 
g cm -3 . We can see that 13 C is depleted at a smaller fraction of 12 C than in Figure [131 
This is because there are more leak nuclei to capture neutrons. The same associated loss of 
accuracy after 13 C depletion is observed. We can also see that the neutron excess evolution 
is closer to 0.3 in this Figure, as expected from equation (151]) . However, as 13 C is depleted, 
the e~-capture rate decreases accordingly and the neutron excess evolution becomes slower. 

Finally, note that although by introducing the dummy leak nuclei we can keep the 
number of independent variables small, we will necessarily lose some accuracy matching the 
neutron excess. Many leak nuclei can capture protons, neutrons or a-particles, but they can 
also capture electrons like 23 Na does. Electron captures on free protons were found to be 
negligible. In order of increasing density, or Fermi electron energy, the following reactions 
may compete with proton, neutron and a-captures: 22 Na(e~, z/ e ) 22 Ne, 24 Na(e~, z/ e ) 24 Ne, 



25 Mg(e", 




(see 


lb en 


1978) 



It is not possible to model the former electron captures accurately without solving 
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for each leak species and its electron capture counterpart independently. Moreover, the 
relative abundances of the different leak species will depend on the temperature and density 
history of the gas, which makes it very difficult to compute the contribution of the different 
electron capture reactions knowing only Yl- However, assuming that the excess of neutrons 
over protons of the leak nuclei is 1.5 under N2, a 10% or 20% accuracy for the neutron 
excess evolution was achieved before and after 13 C depletion, respectively. This translates 
into a maximum of 5% or 10% error for the electron fraction, normally below 1% for small 
quantities of burnt 12 C. 

6.0.1. Timing and Energy Generation Comparison 

In Figure [15] we show the ratio between the elapsed times of the simplified and full 
integrations vs the fraction of burnt 12 C. The Nl and N2 sub-scripts correspond to the 
solution of the first and second simplified networks, respectively. We can see that the Nl 
can over-estimate the speed at which 12 C burns by more than 20%, whereas the second 
simplified network reproduces the time evolution to better than 5% error. This is due to 
the absence of leak reactions in the first approximation, which over-estimates the amount 
of 12 C proton and neutron captures at a given time. 

In Figure [16] we show the ratios between the photon and neutrino energy release rates 
in Nl and N2 and the photon and neutrino energy release rate in the full network. We can 
see that the photon energy release rates is typically off in Nl by 20% or more, whereas 
N2 matches them at the 5% level, except after 13 C is depleted. This is because the leak 
reactions tend to release more energy than the 12 C proton and neutron captures. The 
neutrino energy release rates are not as accurate because secondary electron captures not 
included in these networks can have an important contribution, but since photon energy 
rates are generally much bigger, the error in the net energy release remains at the 5% level. 
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4.0xl0 8 K, 1.0xl0 9 gcm~ 




|AY( 12 C)/Y ( 12 C)| 



Fig. 14. — Same as Figure [TU but changing the density to 10 9 g cm -3 . We can see that 
the network N2 loses its accuracy at a lower fraction of burnt 12 C, which is due to stronger 
neutron leaks which cause 13 C to be relatively depleted sooner. 
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4.0 xlO 8 K, 3.0 xlO 9 gem" 




|AY( 12 C)/Y ( 12 C)| 



Fig. 15. — The ratio between the elapsed time in the two simplified networks and the full 
nuclear network plotted against the fraction of burnt 12 C. Nl and N2 sub-scripts indicate 
which simplified network has been used for the comparison. 
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Fig. 16. — Same as Figure [151 but we now show the ratios of the photon and neutrino energy 
release rates between the simplified and the full networks. Nl and N2 sub-scripts indicate 
which simplified network has been used for the comparison. 
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7. DISCUSSION AND CONCLUSIONS 

We have derived two approximate nuclear networks that can be used in the hydrostatic 
carbon burning regime of carbon-oxygen white dwarfs (CO WDs) approaching ignition. 
The networks have the advantage of being able to track accurately the relevant species 
during this phase of evolution without having to include the fast-evolving protons, neutrons, 
a-particles or 13 N nuclei. 



Using the same integration method (jBader fc Deuflhard 



1983), convergence tolerance 



and initial conditions in one-zone models, we have found that the integration of Nl or N2 is 
much faster than that of the detailed network. Depending on the temperature, density and 
composition, the ratio between the integration times of the detailed network and N2 varied 
between a factor of a few and a factor of several hundreds. Integration time ratios between 
N2 and Nl varied between 1.2 and 1.3. Note that the slowest step in the integration of 
full stellar evolution models is normally Gaussian elimination or similar methods, which 
scale with the cube of the number of independent physical variables. Thus, the simplified 
networks presented in this work could have important applications in full stellar evolution 
models. 

Depending on the amount of 12 C being burnt and the desired errors in the mole 
fractions of the dominant nuclei, the time and the net energy release, either the Nl or N2 
approximations described in Sections |3] and could be used in CO WD interiors. For an 
accuracy of the order of 50%, and to have a simple picture of the main flows involved, Nl 
can be used with caution. For an accuracy of the order of 5%, N2 should be used instead. 
If N2 is used and the criterion in equation (j55"|) does not hold, only an accuracy of the order 
of 20% can be guaranteed, as long as all the trace nuclei time-scales are much shorter than 
the 12 C-burning time-scale. 

These accuracies are not valid for nuclei whose mole fractions are not significant, for 
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example 23 Ne when the /3-decay time-scale is much shorter than the e~-capture time-scale, 
since in this case the 22 Ne(n, 7) 23 Ne reaction will dominate its evolution, a reaction that 
cannot be included accurately under N2. We have also found that for relatively low 
densities, below 10 8 g cm -3 , the quoted errors can double, but only when significant 
amounts of 12 C have been burnt, which normally occurs when the density is significantly 
above this number. Thus, N2 is more accurate for central ignition models, but it could also 
be used in off-center ignition models. 

We have shown how to derive the approximations and have compared them to a detailed 
network at fixed temperature and density. We have also discussed how the approximations 
can break down and when they can be used. Since all the details of the derivation are 
shown, these simplified networks can be further improved straightforwardly. It is also 
possible to build intermediate approximations between Nl and N2, for example, including 
21 Ne in the dummy leak species defined in equation fj4*2|) to remove one independent variable 
from the solution. We have tested this last approximation in a few cases and the resulting 
errors appears to be twice the errors in N2. 

Although for simplicity we have only shown fixed temperature and density integrations, 
these approximations can be used in environments with varying temperature and density, 
such as real stellar evolution models. This is because in the temperature and density range 
found in pre-ignition WD interiors the trace nuclei reach their equilibrium values much 
faster than the typical environmental variables vary inside the star, even within strong 
pre-ignition convective velocity fields. 

The networks can account for 23 Na or 13 N e~-captures for the purposes of following 
the evolution of the pressure-supporting electrons in pre-supernova WDs. We have shown 
that they will be valid even in convective WD interiors. We recommend the use of these 
networks when the mass fractions of the most abundant species are relevant, namely 12 C, 
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16 0, 13 C, 20 Ne, 23 Ne or 23 Na, or the electron mole fraction, Y e , or the energy generation 
rates. We do not recommend its use if the abundances of other nuclei not included in this 
discussion are being studied. 

We have introduced an important tool to understand the effect of the convective Urea 
process on the ignition conditions of SNe la. We foresee the application of these simplified 
networks or their modification in detailed one or multi-dimensional stellar evolution models 
trying to understand pre-supernova CO WDs or si milar objects (see e.g. 



Lesafire et al. 



2006 



Zingale et al. 



2009 



Iapichino fc Lesafire 



2010). 
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